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Abstract. Unwanted stiction in micro- and nanomechanical (NEMS/MEMS) systems due to dispersion 
(van der Waals, or Casimir) forces is a significant hurdle in the fabrication of systems with moving parts on 
these length scales. Introducing a suitably dielectric liquid in the interspace between bodies has previously 
been demonstrated to render dispersion forces repulsive, or even to switch sign as a function of separation. 
Making use of recently available permittivity data calculated by us we show that such a remarkable non- 
monotonic Casimir force, changing from attractive to repulsive as separation increases, can in fact be 
observed in systems where constituent materials are in standard NEMS/MEMS use requiring no special 
or exotic materials. No such nonmonotonic behaviour has been measured to date. We calculate the force 
between a silica sphere and a flat surface of either zinc oxide or hafnia, two materials which are among the 
most prominent for practical microelectrical and microoptical devices. Our results explicate the need for 
highly accurate permittivity functions of the materials involved for frequencies from optical to far-infrared 
frequencies. A careful analysis of the Casimir interaction is presented, and we show how the change in the 
sign of the interaction can be understood as a result of multiple crossings of the dielectric functions of the 
three media involved in a given set-up. 
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1 Introduction 

More than 60 years ago Casimir predicted pQ that bound- 
ary effects on the electromagnetic fluctuations can pro- 
duce attraction between a pair of parallel, closely spaced, 
perfectly conducting plates. His calculation was extended 
to real materials by Lifshitz [2J. Only half a century later 
was precision measurement of the Casimir-Lifshitz force 
between macroscopic bodies made possible [3JE31E51E3 > con- 
firming the theory of Casimir and Lifshitz even though 
certain discrepancies between theory and experiment still 
persist [TJ[8] . Although in all its most immediate manifes- 
tations the Casimir force is attractive, theoretical schemes 
have long existed whereby Casimir repulsion may be achieved 
jMHElGMS]- Munday, Capasso, and Parsegian [H] fa- 
mously demonstrated that the Casimir-Lifshitz force could 
be repulsive by a suitable choice of interacting surfaces in 
a fluid, following similar experiments preceding it |15j . 

We show in the present work that Casimir repulsion, 
and even transitions from attraction to repulsion with 
varying separation, are possible with some of the most 
important materials in use in the fields of micro and na- 
noelectronics and microoptics. Although the introduction 



of a dielectric liquid is still required (we suggest bromoben- 
zene as used in |14j). no further use of special or exotic 
materials is necessary. Change of sign of the Casimir force 
with increasing separation was analyzed by Phan and Viet 
|16) . and by Bostrom et al. [17, 18 and the concept was al- 
ready familiar from theory and experiments in the context 
of films on surfaces (9,19 ( 20,2l|. Here we present calcula- 
tions making use of two important materials within mi- 
croelectronics, zinc oxide (ZnO) and hafnia (HfC^), made 
possible by newly calculated dielectric data for these ma- 
terials. In simplistic terms, the necessary requirement for 
Casimir repulsion is that the permittivities of the two bod- 
ies and the interspatial liquid satisfy the inequality 



£body 1 ^ ^interspace ^ S-body 2- 



(1) 
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Subtleties such as change of sign occurs because the di- 
electric response of a medium to an imposed field depends 
strongly on the field's frequency, so the inequality may be 
satisfied in some frequency ranges, but not in others. The 
Casimir force depends on the response over a broad range 
of frequencies, although the main contributions come from 
frequencies where uid/c < 1 (d is the separation between 
bodies), which explains how the sign of the force can de- 
pend on the body-body separation. 
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We have calculated the Casimir force between materi- 
als which are already important in nano- and microelec- 
tromechanical systems (NEMS and MEMS), separated by 
the oil bromobenzene, to show that these standard materi- 
als are sufficient to observe both attraction and repulsion 
in one and the same set-up. Zinc oxide is a multi-functional 
semiconductor much used in optoelectronic devices. The 
material is transparent at optical frequencies while block- 
ing the ultraviolet light, and it becomes highly conductive 
when n-type doped. Hafnia is commonly used for opti- 
cal coating and is a leading candidate for the replacement 
of silica for a number of microelectronic applications due 
to its high permittivity, low optical absorption, and low 
thermal expansion |22j . 

To make contact with experiments [14] we consider 
the geometry consisting of a silica sphere (such as may 
be attached to an atomic force microscope tip [3]) inter- 
acting with a ZnO or hafnia surface across bromobenzene 
(Bb). Apart from demonstrating the possibility of transi- 
tion from attractive to repulsive Casimir forces, our cal- 
culations demonstrate the importance of having access to 
high accuracy dielectric functions [23*11241125] . Different lev- 
els of modeling of the dielectric functions give radically 
different results. For instance, the sign of the force be- 
yond about 50 nanometers may even change depending 
on whether low frequency electron-phonon contributions 
are included in the modeling or not. We review the cal- 
culations of the dielectric functions of the different me- 
dia and recapitulate the theory of Casimir-Lifshitz forces, 
whereupon we present numerical results. 



2 Calculation of dielectric permittivities 
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Fig. 1. (Colour on-line) Dielectric functions at imaginary fre- 
quencies for ZnO, Si02, Bb (bromobenzene), and Hf02. The 
results for ZnO and Hf02 are from calculations in this work us- 
ing the GW method based on the DFT [281132] : the first ZnO 
version neglects and the second includes phonon contributions. 
For Si02 the first version is from an oscillator model [26], while 
the second is the present GWo results. For Bb the first result 
[14] . is from an oscillator model; the second [33] from a differ- 
ent oscillator model; note that for the second of the dielectric 
functions for Bb the static value only entering the n = term 
is 5.37. 



The Casimir force at temperature T can be calculated if 
the dielectric functions (for discrete imaginary frequencies, 

= 2irnk&T /h) are known. The dielectric functions of 
the materials play the essential role in the Lifshitz theory 
[9] . Figure [T] presents the dielectric functions on the imag- 
inary frequency axis for bromobenzene (Bb) determined 
in [T3], Si02 (both calculated and a modeled dispersion 
in [5S] based on experiments), ZnO, and Hf02- For the 
theoretically determined dielectric functions in this work, 
we present the calculations both including and exclud- 
ing the optical phonon modes in order to illustrate the 
importance of modeling the electron-phonon coupling for 
analyzing the spectra at frequencies below ~ 10 15 rad/s. 

The complex dielectric functions for wurtzite ZnO, mon- 
oclinic HTO2 and a-quartz Si02 were determined employ- 
ing a first-principles approach within the density func- 
tional theory (DFT). The electronic structure, neglecting 
electron-phonon coupling, was obtained from the partial 
self-consistent GWq method where the Green functions are 
updated iteratively whereas the screened Coulomb poten- 
tial W is fixed [271125], The imaginary part of the dielec- 
tric function was calculated in the long wavelength limit 
from the joint density of states and the optical momen- 
tum matrix. Since the dielectric function in polar materi- 
als can depend strongly on the electron-phonon coupling, 



we modeled this contribution to the dielectric function us- 
ing the Lorentz model and Kramers-Heisenberg formula 
with multi-phonon contribution [29,30 . The phonon part 
of the dielectric function is: 

^-}£phonon(W) — > J o n : ■ (2 J 

Here, (jJloj and wto.j are the longitudinal optical (LO) 
and transverse optical (TO) phonon frequencies of the jth 
mode, respectively, and Eooj is the high frequency dielec- 
tric constant of the jth phonon mode. We determine eoo,j 
by employing experimental data of the phonon frequen- 
cies [31]. The low-energy spectra is verified by calculat- 
ing the static dielectric constant from the Born effective 
charges. The calculated dielectric functions of ZnO, Si02, 
and HTO2 on the imaginary frequency axis are shown in 
Fig. [T] in the small damping limit (i. e., phonon damping 
parameters 7^ — » ). The corresponding static dielectric 
constants are for ZnO 7.9 (8.1-8.3), for Si0 2 3.9 (3.9-4.4), 
and for Hf02 24 (15-25), respectively. The numbers within 
parentheses are the experimental values from Refs. [34,35, 

05II37J- 
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3 The Casimir force formalism 

Using the Deryaguin (or proximity force) approximation 
[3"8] the Casimir-Lifshitz force of a planar surface of mate- 
rial 1 (ZnO or HTO2) with a sphere (radius R) of material 
3 (S1O2) across medium 2 (Bb) results in a summation of 
imaginary frequency terms [9"1 I31?II3"8] : 



F = 2nR 



OO 

n=0 



(3) 



Note that positive values of F correspond to repulsive 
force. 

In the retarded treatment there are contributions from 
the two light polarisations, transverse magnetic (TM) and 
transverse electric (TE), g(£ n ) = g™(£„)+g TE (£„), where 



--k B T 
x In 

~-k B T 
x In 



d 2 q 
(2tt)* 
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d 2 q 
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and the Fresncl reflection coefficients are 

£a(^n)7i(^n) - Ei(ifn)72(«£n) 
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. (4) 

(5) 
(6) 

(7a) 
(7b) 



The non-retarded limit can be investigated by letting the 
velocity of light go to infinitjQ: 



F,. 



Rk B T 
4d 2 



n=0 



'u 3 [r™(iZ n )r^(itn)], (8) 



where the order 3 polylogarithm is 
Li 3 (x) 



Jl. 

n=l 



(9) 



The non-retarded approximation is the short separation 
asymptote of the real (retarded) Casimir force. 

Frequency intervals where the intervening medium has 
a dielectric permittivity in between the permittivities of 
the two bodies give a repulsive contribution as t^it^ < 
for both polarisations in these intervals; other intervals 
give an attractive contribution. Correct calculation of the 
force — even its sign — thus requires accurate dielec- 
tric functions for a wide frequency region. Since, roughly 



In the model case where reflectivity does not tend to zero 
as £ — 5- 00, it is necessary to retain finite speed of light for 
convergence. For perfectly reflecting surfaces, thus, there is no 
non-retarded regime. 





in 




LO ! 


Is 










11) 






« 




s 


11) 




£r 






10 ~ s 




10 ~' 




10- 1 




-10~ 


■g/cm 


-UT 


to 




% 


-10~ 












-UT 




-111 



Retarded force 



; Retarded force 



Only n — term 



Non-retarded force 



Separation (A) 

Fig. 2. (Colour on-line) Casimir-Lifshitz force between ZnO 
(1) surface and Si02 (1) sphere in Bb (1). The fully retarded 
free energy is attractive in the small and large separation limits, 
but is repulsive in an intermediate interval. The nonretarded 
energy and the n = contribution are both attractive and fall 
off with increasing d as d~ 2 . 



speaking, the expression ((3J) picks up its main contribution 
from £ < c/d, the sign of the force in the short and long 
separation regimes can be designed by choosing materi- 
als with the appropriate dielectric properties in different 
frequency ranges. In particular, the long-separation range 
d ;§> c/£i is dictated entirely by the n = term of the 
sum, 



_ Rk B T 
F n=0 - — gjg-Lia 



e 2 (0)-e 3 (0) e 2 (0)-ei(0) 



£ 2 (0)+e 3 (0) e 2 (0)+ei(0) 



(10) 



Since we are limited by the proximity force approximation 
to d <C R this long-separation range is only valid if R is 
large enough. In Fig. [2] we see that for our choice of mate- 
rials the n = asymptote starts to dominate the result at 
fractions of a micrometer. Thus, for the long-separation 
range to be of interest here the radius of the sphere has to 
be of micrometer size or larger. To avoid misunderstand- 
ings we repeat that all presented results are valid only for 
d <C R. For larger separations, d > R or larger, all results 
are different. The n — term is different and represents 
the long-separation asymptote that is approached when 
d — y oo. 



4 Numerical examples 

Using the theory and permittivity data laid out in the 
above sections, numerical calculations are straightforward. 
All results presented here are for 300 K. We show in Figs. [5] 
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Fig. 3. (Colour on-line) Casimir-Lifshitz force between ZnO 
(2) surface and SiCh (2) sphere in Bb (2). The nonretarded 
energy is attractive whereas the n = contribution is repulsive, 
whereas the fully retarded (actual) force changes sign at an 
intermediate separation. Both asymptotes decrease as d~ 2 . 
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Fig. 4. (Colour on-line) Casimir-Lifshitz force between Hf02 
surface and Si02 (2) sphere in Bb (2). The nonretarded energy 
is attractive while the n = contribution is repulsive, and the 
fully retarded (actual) force tends to these asymptotically in 
the short and long separation limits, respectively. 



and[3]the retarded force, nonretarded force, and the n = 
term. It is clear that retardation changes in a major way 
the Casimir-Lifshitz force between a silica sphere and zinc 
oxide surface in bromobenzene. The figures are in log-log 
scale so as to visualise the changing sign of the forces. 

By means of optical measurements or accurate theo- 
retical calculations of the dielectric properties of the in- 
teracting objects and liquid it is possible to predict the 
force — from short-range attractive van der Waals force 
to intermediate range repulsive Casimir force. The long 
range entropic asymptote is either attractive or repulsive 
depending on the optical properties for low frequencies. 

Upon inspection of figures [5] and [3] it is striking that 
the two different permittivity models can predict different 
signs for the large separation asymptote of the Casimir 
force between silica and ZnO. Mathematically this can be 
understood from equation (ITU1) by regarding the zero fre- 
quency limit of the different permittivities involved: the 
low frequency asymptote of the permittivity of the liquid 
bromobenzene falls either between or above those of the 
solids depending on whether or not low frequency electron- 
phonon contributions are included for ZnO. This high- 
lights that in order to predict the Casimir force at larger 
separations, indeed even its sign, care must be taken that 
the model permittivities used take realistic values in the 
quasistatic limit. Note that the predictions in the two cases 
differ greatly at all separations above about 100 nm, where 
the Casimir effect can still play a practical role in MEMS, 
thus the drastic effects of having inadequate permittiv- 
ity data is not just a curiosity but could have important 
consequences. 

Figures [3] and 0] present our best estimates for the 
Casimir effect between a silica sphere and, respectively, 
ZnO and Hf02, in bromobenzene. These calculations draw 
on the best available dielectric data for ZnO, Si02 and 
Hf02 as determined from the GWo calculations, whereas 
permittivity data for Bb are taken from [33 . Both of 
these systems, which apart from the interspatial oil con- 
sist only of typical MEMS materials, exhibit similar be- 
haviour: the Casimir force is attractive at short separa- 
tions and becomes repulsive at larger separations. No such 
non-monotonic Casimir force has been measured to date, 
yet our calculations predict that it may be not only ob- 
servable, but even of practical importance in realistic mi- 
croelectromechanical set-ups. The proximity force approx- 
imation is valid if the sphere radius is much larger than 
the separation [40ll41j . Thus the sphere radius must be in 
the micro meter range or larger for the rightmost parts of 
figures f2]|3] to be valid. 



5 Conclusions 

We have shown herein that non-monotonic Casimir force, 
changing from attractive to repulsive with increasing sep- 
aration, may be found in set-ups in which the constituent 
materials are in standard use in NEMS and MEMS. To 
wit, we have calculated the Casimir force between a silica 
sphere and a half-space of ZnO and Hf02, respectively, 



M. Bostrom et al.: Casimir attractive-repulsive transition in MEMS 



5 



immersed in bromobenzene. All three solids are among 
the primary materials used in NEMS and MEMS. 

The striking change of sign of the Casimir force is a 
consequence of the frequency dependence of the three dif- 
ferent permittivities involved. At small separations, all fre- 
quencies contribute to the force, whereas the force at long- 
separations depends on the materials' quasistatic dielec- 
tric response only. It is possible therefore for the force to 
change sign more than once, and the force at long sep- 
arations, indeed even its sign, depends sensitively on the 
low frequency asymptotics of the dielectric function model 
used. It is of vital importance for accurate force calcula- 
tion that low frequency contributions from the far-infrared 
regime, such as photon-phonon interactions, are included 
in the permittivity function. This frequency regime is typ- 
ically not covered in tables of optical data, but its exclu- 
sion could cause radically wrong predictions of the Casimir 
force, not only in the long-separation asymptote, but also 
at transitional separations between the short and long dis- 
tance limits. 

We suggest, in conclusion, that measurements of the 
attraction-repulsion transition of the Casimir force may 
not only be possible, but of considerable practical interest. 
In our example systems, seen in figures [3] and [4j the transi- 
tion happens at about 150 nm. This transition separation 
can be modified by slighly altering the dielectric proper- 
ties involved. It is straightforward to generalize the results, 
for example, by using ultra-thin coatings of, e.g., molyb- 
denum disulfide or graphene [4"21l4"3"] . Such surface modi- 
fications may alter the way Casimir-Lifshitz forces switch 
from attraction to repulsion to attraction/repulsion. 

M.B. acknowledges support from an ESF exchange grant within 
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fect" , through the network CASIMIR. CP. and M. B. acknowl- 
edge support from VR (Contract No. 90499401) and STEM 
(Contract No. 34138-1). B.E.S. acknowledges financial support 
from VR (Contract No. 70529001). 
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